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ABSTRACT 


In  this  paper  we  address  a  basic  problem  In  machine  perception:  the 
tracing  of  "llne-llke”  structures  appearing  In  an  Image.  It  Is  shown 
that  this  problem  can  profitably  be  viewed  as  the  process  of  finding 
skeletons  In  a  gray  scale  Image  after  observing  (1)  that  line  detection 
does  not  necessarily  depend  on  gradient  Information,  but  rather  Is 
approachable  from  the  standpoint  of  measuring  total  Intensity  variation, 
and  (2)  that  smoothing  the  original  Image  produces  an  approximate 
distance  transform.  An  effective  technique  for  extracting  the 
delineating  skeletons  from  an  Image  Is  presented,  and  examples  of  this 
approach  using  aerial,  Industrial,  and  radiographic  Imagery  are  shown. 


I  INTRODUCTION 


For  many  tasks  In  scene  analysis,  there  may  not  exist  general 
solutions  Independent  of  purpose  or  Intended  application.  However,  for 
the  task  of  linear  delineation,  one  can  easily  find  Image  subsets  for 
which  a  panel  of  human  observers  would  be  almost  unanimous  In  their 
Interpretation  without  having  to  agree  on  the  explicit  criteria 
underlying  their  decisions;  our  goal  Is  to  produce  a  computer  system 
that  can  perform  the  delineation  task  at  close  to  human  levels  for  at 
least  these  more  obvious  cases,  especially  where  semantic  knowledge  Is 
not  required.  In  this  paper  we  present  some  new  ways  of  looking  at  the 
problem  of  linear  delineation  and  provide  techniques  that  are 
significantly  more  general  and  effective  than  previously  reported 
methods  for  this  task. 


II  PROBLEM  DEFINITION 


For  the  purposes  of  this  paper,  we  define  linear  delineation  as  the 
task  of  generating  a  set  of  lists  of  points,  for  a  given  2-D  Image,  such 
that  the  points  In  each  list  fall  sequentially  along  what  any  reasonable 
human  observer  would  describe  as  a  clearly  visible  "llne-llke"  structure 
In  the  Image.  Practical  examples  of  this  task  might  be  to  delineate  the 
roads,  rivers,  and  rail  lines  In  an  aerial  photograph,  or  to  trace  the 
paths  taken  by  blood  vessels  In  a  radiographic  angiogram,  or  to  locate 
the  wiring  paths  on  a  printed  circuit  board;  however,  our  goal  In  this 
paper  is  not  to  look  for  specific  real-world  objects  or  to  assign 
semantic  labels  to  the  detected  linear  structures,  but  rather  to  find 
what  a  human  observer  might  choose  as  the  most  perceptually  obvious 
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occurrences  of  such  structures.  We  further  distinguish  between  the 
problems  of  (1)  detecting  the  edges  or  contours  of  extended  objects,  and 
(2)  delineating  those  objects  whose  appearance  is  adeqiiately  represented 
by  a  central  skeleton  —  only  the  second  problem  is  addressed. 


Ill  LINES  AND  EDGES 


While  most  approaches  to  linear  delineation  do  not  distinguish 
between  lines  and  edges,  and  even  use  edge  detection  as  a  necessary 
first  step  in  the  delineation  task,  a  critical  concept  advanced  in  this 
paper  is  the  distinction  between  line  and  edge  detection. 

Edge  detection  is  based  on  the  concept  of  finding  a  discontinuity 
(in  intensity  or  some  other  locally  measurable  attribute  such  as  color 
or  texture)  between  two  adjacent  but  distinct  regions  in  an  image. 
However,  in  a  digital  representation  of  an  image,  a  smooth  surface  can 
always  be  fit  to  the  sample  values  of  the  Integer  raster.  Thus,  edge 
detection  must  be  based  on  parameters  or  thresholds  set  by  assumptions 
about  the  nature  of  the  image.  Even  if  edge  points  are  marked  only  at 
those  locations  at  which  there  are  first  derivative  maxima,  or  zeros  of 
the  second  derivative,  ultimately,  an  arbitrary  decision  is  made  in 
deciding  when  the  corresponding  gradient  is  large  enough  to  be  called  a 
discontinuity. 

Intuitively,  a  "llne-llke"  or  linear  structure  is  a  (connected) 
region  that  is  very  long  relative  to  its  width,  and  has  a  ridge  or 
skeleton  along  which  the  intensities  change  slowly  and  are  distinguished 
from  intensities  outside  the  region;  the  width  need  not  be  constant,  but 
any  changes  in  width  should  occur  in  a  smooth  manner.  To  simplify  the 
discussion,  we  will  assume  the  linear  structures  are  distinguished  by 
their  rldgepolnts  being  brighter  than  the  surrounding  background,  but 
any  other  specified  attribute,  which  is  locally  detectable,  would  be  an 
acceptable  substitute.  It  is  Important  to  recognize  the  fact  that  a 
clearly  visible  line  in  an  image  may  not  have  locally  detectable  edges 
and  thus  no  locally  measurable  width,  or  possibly  only  one  detectable 
edge,  or  even  two  edges  which  are  significantly  separated  and 
nonparallel  as  might  occur  in  a  local  widening  of  a  river.  It  is  also 
generally  the  case  that  linear  structures  have  no  visible  Internal 
detail  that  is  essential  to  their  delineation. 

As  a  point  of  interest,  it  might  be  noted  that  the  mechanisms  for 
generating  subjective  edge  and  line  illusions  are  quite  different; 
subjective  edges  appear  to  require  a  3-D  interpretation,  while 
subjective  lines  appear  to  be  produced  by  adaptation  phenomena. 
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IV  SMOOTHING,  DISTANCE  TRANSFORMS,  AND  THE  GRAY  SCALE  SKELETON 


If  we  can  find  the  edges  of  a  linear  structure,  we  can  generate  a 
distance  transform  and  extract  a  skeleton  as  the  desired  delineation 
(e.g.,  Rosenfeld  [1],  Flschler  {2]).  However,  as  noted  In  the  preceding 
section,  linear  structures  do  not  always  have  locally  detectable  edges, 
and,  since  all  of  the  generally  known  techniques  for  deriving  a  skeleton 
require  a  complete  contour,  some  other  approach  Is  required.  The 
classical  skeletonizing  techniques  Intimately  link  the  contour/edges  of 
a  region  and  Its  skeleton,  and  It  Is  Just  this  linkage  that  we  wish  to 
break. 

Surprisingly,  something  equivalent  to  a  distance  transform  that 
works  on  gray  scale  Images,  and  on  binary  Images  as  well.  Is  already 
available.  To  achieve  our  purpose,  we  need  only  observe  that  the 
Intensities  In  a  properly  smoothed  Image  can  be  considered  to  be  the 
values  of  an  approximate  distance  transform.  What  Is  the  best  smoothing 
function  for  general  use?  Actually,  It  doesn't  seem  to  make  much 
difference  In  many  cases.  Most  digital  Images  have  been  processed  by 
low-pass  optical  and  electronic  systems  that  have  Inserted  the  required 
minimum  level  of  smoothing.  The  viewpoint  that  the  smoothed  Image  can 
be  considered  to  be  a  distance  transform  Is  the  essential  element. 
However,  If  we  start  with  a  binary  Image,  or  a  very  noisy  Image,  then 
additional  smoothing  la  often  desirable.  Since  we  are  not  concerned 
with  blurring  edges,  and  we  would  like  to  eliminate  or  blur  any 
structure  or  texture  Internal  to  the  linear  regions,  we  want  the 
smoothing  function  to  have  a  width  of  approximately  that  of  the  region 
to  be  delineated.  If  the  width  of  the  smoothing  function  Is  Increased 
further,  the  thinner  linear  structures  are  eventually  eliminated.  Thus, 
If  we  wish  to  find  all  possible  linear  structures  without  prior 
knowledge  of  the  content  of  the  Image,  the  processing  should  be  repeated 
with  a  set  of  smoothing  filters  having  a  spectrum  of  widths.  Actually, 
no  more  than  two  or  three  filtering  steps  should  ever  be  required.  For 
example,  to  trace  all  the  linear  structures  (diameters  up  to  20  pixels) 
In  a  noisy  radiographic  angiogram,  a  single  filter  of  width  20  was  used 
(see  Figure  1).  Smoothing  Introduced  by  the  acquisition  process  was 
sufficient  to  produce  excellent  results  In  tracing  the  linear  structures 
In  aerial  Imagery  (see  Figures  4  and  5). 


V  RIDGES  (OR  VALLEYS),  OPERATORS,  AND  NEIGHBORHOODS 


Having  produced  an  approximate  distance  transform  via  smoothing,  we 
now  must  deal  with  the  problem  of  locating  the  rldgepolnts  that  denote 
the  spines  (skeletons)  of  the  linear  structures.  When  an  exact  distance 
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transform  is  derived  from  a  complete  contour,  noise  is  not  a  problem  and 
the  skeleton  has  assured  geometric  properties  that  make  it  easy  to 
detect;  finding  the  rldgepoints  of  an  approximate  distance  transform  is 
considerably  more  complex. 

We  traditionally  distinguish  between  locally  and  globally 
detectable  features:  local  features  are  detectable  by  an  intensity 
pattern  which  can  be  observed  through  a  small  peep-hole  centered  on  the 
feature,  while  global  features  are  ambiguous  in  a  small  area.  The  model 
or  description  of  the  local  feature  is  generally  compiled  into  an 
intensity  patch  (matched  filter  or  operator)  which  can  be  convolved  with 
the  image  to  detect  the  corresponding  feature.  In  the  case  of  an  exact 
distance  transform,  a  3X3  pixel  operator  is  sufficient  to  detect 
rldgepoints  (a  2X2  operator  is  sufficient  for  the  Labeled  Distance 
Transform  (Flschler  [2]));  for  the  approximate  distance  transform,  a 
small  fixed-size  operator  is  ineffective. 

The  principal  utility  of  a  local  operator  is  that  the  number  of 
data  patterns  the  operator  might  encounter  is  small  enough  to  allow  one 
to  enumerate  a  decision  for  each  such  pattern.  If  we  further  agree  to 
use  a  small  square  window  of  the  image  as  our  local  domain,  and  to  use 
either  table  look-up  or  convolution  as  the  basis  for  decision  making, 
then  a  uniform  mechanization  can  be  employed  to  Implement  a  large  number 
of  distinct  (and  generally  unrelated)  local  operators.  The 
attractiveness  of  this  second  Implementation  aspect  has  led  to  the 
situation  that  almost  all  low  level  (local)  scene  analysis  is  done  using 
such  peep-hole  type  operators.  The  disadvantage  of  this  approach  is 
that  the  concept  of  local  is  relative  to  the  size  of  the  entity  of 
Interest,  and  either  one  must  know  this  size  in  advance,  or  use  a  whole 
family  of  operators  of  increasing  size,  where  the  larger  operators  lose 
the  advantages  that  led  to  their  use  in  the  first  place.  In  the  case  of 
line  detection,  where  the  line  width  can  vary  over  a  wide  range  of 
values,  the  conventional  operator  concept  is  inappropriate. 

Based  on  these  general  issues  (even  more  than  on  the  Immediate 
problem  at  hand),  we  have  considered  other  realizations  of  general 
"local”  decision-making  processes  that  satisfy  the  previously  stated 
conditions,  but  do  not  necessarily  lend  themselves  to  a  convolution  type 
mechanization.  In  particular,  restricting  our  attention  to  finding  the 
maxima  and  minima  of  functions  of  the  displacement  along  a  space  curve 
defined  over  the  image,  satisfies  our  requirements  for  computational  and 
decision-making  simplicity  even  when  the  curve  traverses  the  entire 
image.  While  the  space  curve  might  assume  any  shape  (e.g.,  follow  the 
contour  of  an  object),  the  analysis  itself  is  independent  of  the  shape; 
for  the  linear  delineation  problem,  we  used  image  intensity  as  a 
function  of  displacement  along  horizontal  and  vertical  scan  lines. 
Since  maxima  and  minima  are  symmetrical  attributes,  we  will  only  discuss 
the  problem  of  labeling  maximal  points  along  the  curve. 
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The  problem  of  finding  the  rldgepolnts  of  an  approximate  distance 
transform  can  be  viewed  as  the  problem  of  finding  the  rldgepolnts  (local 
maxima)  of  an  exact  distance  transform  to  which  some  amount  of  noise  has 
been  added.  We  are  not  concerned  about  the  possibility  of  making 
Isolated  (Incoherent)  Incorrect  decisions,  because  we  have  developed 
effective  linking  and  pruning  methods,  described  In  the  next  section, 
that  are  capable  of  eliminating  such  errors.  Our  main  problem  Is  that 
we  cannot  count  on  finding  either  large  local  gradients  or  using  known 
line  width  to  determine  some  minimum  significant  gradient  threshold  to 
Identify  valid  rldgepolnts;  additionally,  noise  will  Introduce  many 
false  local  maxima.  Thus,  we  must  use  total  Intensity  change,  rather 
than  rate  of  change,  to  detect  valid  rldgepolnts,  and  we  must  have  an 
effective  way  of  determining  such  total  change  even  In  the  presence  of 
local  variation  Introduced  by  noise.  (While  It  Is  not  Immediately 
obvious  that  total  Intensity  change,  rather  than  rate  of  change,  will 
recover  the  perceptually  obvious  linear  features,  our  experiments 

indicate  that  this  Is  Indeed  the  case.) 

*  \ 

Our  approach  is  to  evaluate  two  attributes  of  each  of  the  detected 
Intensity  maxima  along  the  space  curves  (in  this  case,  horizontal  and 
vertical  scan  lines),  which  we  call  the  "local"  and  "global"  maxima 
values.  The  local  maxima  value  of  a  point  is  the  total  Intensity 
difference  from  the  point  to  the  highest  of  its  immediate  left  and  right 
intensity  minima  along  the  curve.  The  left  (right)  global  maxima  value 
of  a  point  Is  the  total  intensity  difference  from  the  point  to  the 
lowest  intensity  value  found  moving  to  the  left  (right)  prior  to 
encountering  a  point  with  an  intensity  value  greater  than  that  of  the 
given  point;  the  global  maxima  value  of  the  point  is  the  smaller  of  Its 
left  and  right  global  values.  In  the  case  of  a  plateau,  only  the  center 
point  Is  treated  as  a  maximal  point  and  evaluated  as  previously 
described.  If  a  point  (or  points)  on  a  plateau  has  an  Immediate 
neighbor  with  a  higher  Intensity  value,  it  Is  not  a  maximal  point  and  It 
Is  not  assigned  either  a  local  or  global  value  (actually,  for 
Implementation  purposes,  non-maxlmal  points  are  assigned  zero  values); 
on  the  other  hand,  every  maximal  point  will  have  both  a  local  and  global 
value  where  the  global  value  equals  or  exceeds  the  local  value. 
Figure  2  provides  some  examples  illustrating  the  operations  just 
described. 

We  have  been  proceeding  under  the  assumption  that  a  large  local 
intensity  maxima  (LIM)  denotes  a  significant  event,  but.  In  the  presence 
of  large  variations  in  image  intensity  or  noise,  the  global  Intensity 
maxima  (GIM)  would  be  a  better  detector  of  significant  intensity 
variation;  however.  In  a  well-smoothed  or  relatively  noise-free  image, 
there  might  be  very  little  difference  in  the  Information  contained  in 
the  LIM  and  GIM  measures.  There  Is  also  the  Issue  of  deciding  what  Is  a 
large-enough  value,  of  either  the  LIM  or  GIM,  to  Indicate  significance. 
In  our  unsmoothed  Image  data,  about  1/3  of  the  points  were  maximal 
points,  and,  In  a  smoothed  Image,  this  percentage  Is  much  smaller  (see 
Table  1).  Given  the  linking  and  pruning  techniques  we  describe  in  the 
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next  section.  It  might  be  possible  to  return  all  the  maximal  points  In  a 
binary  mask  and  still  extract  the  desired  line  structure  from  the 
background  noise  contained  In  such  a  mask.  However,  It  makes  much  more 
sense  to  first  eliminate  those  maximal  points  that  do  not  have  enough 
Intensity  variation  to  be  perceptually  distinguishable  from  a  flat 
background.  (It  would  even  appear  that  we  could,  without  losing 
essential  Information,  eliminate  those  maximal  points  with  a  total 
variation  less  than  that  required  to  perceive  them  against  a  random- 
noise  field  with  the  same  statistical  variation  as  the  measured 
variation  over  some  surrounding  neighborhood  In  the  Image.) 

Rather  than  attempting  to  find  some  optimal  threshold  setting  (In 
the  sense  of  maximum  noise  elimination  without  any  loss  of  linear 
structure),  which  would  be  difficult  or  Impossible  to  automatically 
determine  at  this  level  of  Information  organization,  we  Iteratively 
adjust  our  threshold  settings  to  satisfy  a  constraint  based  on  a 
complexity  measure.  These  program-determined  thresholds  typically  allow 
at  least  two  to  ten  times  the  number  of  rldgepolnts  (maxima)  to  be 
retained  above  that  which  would  result  by  manually  setting  the 
thresholds  to  achieve  visually  acceptable  results.  The  final 
elimination  of  "non-slgnlflcant"  maxima  Is  achieved  later  In  the 
processing  at  a  higher  level  of  organization. 


VI  CLUSTERING,  LINKING,  PRUNING,  NODE  ANALYSIS,  RANKING,  AND  FINAL 
DELINEATION 


Based  on  the  availability  of  a  binary  overlay  depicting  the 
locations  of  the  major  linear  structures  contained  In  the  given  gray 
scale  image,  obtained  as  described  In  the  preceding  section,  we  have 
been  able  to  demonstrate  that. 

(1)  The  linking  step  In  the  delineation  process  can 
effectively  be  based  on  the  single  attribute  of  geometric 
proximity,  and  that  a  clustering  or  association  step, 
followed  by  the  construction  of  a  Minimum  Spanning  Tree 
(MST)  through  the  points  of  each  cluster  (Flschler  [3], 

Zahn  [4]),  will  correctly  link  the  rldgepolnts  along  the 
skeletons  of  the  linear  structures. 

(2)  The  desired  delineations  will  be  embedded  In  trees 
containing  additional  branches  that  are  either  minor 
linear  structures  or  noise,  and  that  simple  pruning 
techniques  can  eliminate  most  of  this  unwanted  detail 
(see  Figure  3;  note  that  tree  pruning  can  effectively 
achieve  simplifications  that  would  be  difficult,  If  not 
Impossible,  at  lower  levels  of  organization  of  the 
Information) . 
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(3)  Having  properly  linked  the  rldgepolnts  and  pruned  some  of 

the  smaller  branches  of  the  resulting  trees,  we  can 
extract  long  coherent  paths  by  a  decision  procedure 
applied  at  each  node  of  each  tree.  This  decision 

procedure,  based  on  the  local  branch  attributes  of 
Intensity,  connectivity,  and  directionality  at  each  node, 
assigns  path  connectivity  through  a  node  by  splitting  off 
Incompatible  branches;  any  remaining  ambiguities  (more 
than  two  branches  entering  a  node)  are  resolved  by 
choosing  those  pairings  that  result  In  the  longest  paths. 

(4)  The  paths  obtained  In  the  tree  partitioning  step  can  be 
rank  ordered  with  respect  to  perceptual  quality  by  a 
metric  based  on  the  path  attributes  of  total  length, 
contrast,  and  continuity. 

Details  of  the  procedures  discussed  In  this  section  are  presented 
In  the  appendix. 


VII  EXPERIMENTAL  RESULTS 


Figures  4-10  present  the  results  of  applying  the  delineation 
algorithm  to  aerial.  Industrial,  and  medical  Images.  Our  goal  was  to  be 
able  to  take  Imagery  from  arbitrary  domains,  and  without  any  human 
Intervention  (e.g.,  parameter  adjustment  or  attention  focusing),  produce 
high  quality  delineation  of  the  obvious  linear  structures.  The  results 
show  that  we  have  accomplished  much  of  what  we  originally  Intended.  The 
delineations  achieved  by  the  uniform  parameter  settings  are  quite  good 
In  all  the  Images  with  the  exception  of  the  angiogram  (which  Is 
extremely  noisy,  and  does  not  really  satisfy  the  criteria  of  having 
’’clearly  visible”  linear  structure);  and  even  here,  by  making  an 
appropriate  selection  of  two  parameters,  the  smoothing  diameter  and  the 
association  distance  (required  In  the  clustering  step),  we  obtain  very 
good  results.  By  using  the  values  produced  by  the  ranking  step  of  the 
delineation  algorithm,  we  believe  It  should  be  possible  to  automatically 
search  the  parameter  space  (of  approximately  10  to  100  parameter 
combinations)  to  optimize  the  processing  for  any  given  Image.  We  are 
currently  Investigating  this  possibility. 
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VIII  CONCLUDING  COMMENTS 


We  have  presented  the  viewpoint  that  the  problem  of  delineating  the 
obvious  linear  structures  in  an  image  is  distinct  from  that  of  finding 
edges  or  contours,  and  is  best  viewed  as  the  process  of  finding 
skeletons  in  gray  scale  Images  (i.e.,  that  line  detection  does  not 
necessarily  depend  on  gradient  information,  but  rather  is  approachable 
from  the  standpoint  of  detecting  total  intensity  variation);  as  a 
necessary  step  in  this  process,  we  have  suggested  that  an  approximate 
gray  scale  distance  transform  can  be  attained  by  smoothing  the  original 
image.  We  have  described  an  effective  technique  for  finding  rldgepolnts 
(points  on  the  delineating  skeleton),  and,  in  the  process,  raised  some 
Important  questions  about  the  conventional  approach  to  designing  "local 
operators . " 

The  competence  of  the  rldgepolnt  algorithm  to  abstract  the  linear 
structure  of  an  Image  is  apparent  by  inspection  of  Figures  4  through  10, 
This  effectiveness  is,  in  a  sense,  a  "psychological  discovery."  It 
would  appear  that  rldgepolnts  are  important  perceptual  primitives  which 
may  play  a  significant  role  in  a  variety  of  other  tasks  (e.g., 
perception  of  surface  shape).  A  significant  feature  of  "rldgepolnts"  is 
that  their  locations  are  Independent  of  any  monotonlc  intensity 
transformation  of  the  image,  and  their  geometric  configuration  is 
Independent  of  a  change  of  scale.  These  are  precisely  the  properties  we 
would  expect  of  a  perceptual  primitive,  but  are  lacking  in  the  commonly 
employed  gradient  dependent  primitives,  such  as  "edgepoints . " 

Starting  with  both  the  binary  overlay  (produced  as  discussed  in  the 
main  body  of  this  paper)  and  the  original  gray  scale  image,  we  have 
demonstrated  via  examples  that  the  remaining  steps  in  the  delineation 
process  can  be  effectively  achieved. 

Our  goal  in  this  work  has  been  to  approach  human  levels  of 
performance  in  finding  perceptually  obvious  delineations  in  Images 
selected  at  random  from  a  reasonably  broad  class  of  scene  domains,  and 
without  any  human  intervention  or  prior  knowledge  about  the  image 
content.  We  believe  that  this  goal  can  be  achieved  through  extension 
and  refinement  of  the  techniques  described  in  this  paper. 
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Appendix  A 


This  appendix  describes  the  procedures  Invoked  by  the  delineation 
algorithm.  These  Include  : 

a)  Smoothing 

b)  Rldgepolnt  Detection 

c)  Association  (grouping  or  clustering)  and  Thresholding 

d )  LI nklng 

e)  Pruning 

f)  Linear  Segment  Extraction 

g)  Linear  Segment  Decomposition 

h)  Linear  Segment  Evaluation 


SMOOTHING 

a)  Purpose:  enhance  connectivity,  especially  of  "thick"  lines,  In 
the  presence  of  noise  (e.g.,  along  noisy  plateaus  —  see  Figure  11). 

b)  Algorithm:  2-D  Gaussian  smoothing  approximated  by  a  1-D 
binomial  mask  Iteratively  convolved  with  the  Image  In  the  horizontal  and 
vertical  directions  (e.g.,  a  mask  of  diameter  five  Is  the  sequence 
<1464  1>). 

c)  Parameter  Space:  mask  diameter  Is  the  only  adjustable 
parameter.  Values  of  <0  3  5  10  20>  would  probably  be  adequate  for  any 
256X256  Image.  In  our  experiments,  we  used  values  of  0  or  3  for  all  the 
Images  except  the  very  noisy  angiogram  which  required  a  smoothing 
diameter  of  20  for  best  results. 

d)  Parameter  Selection:  since  smoothing  can  alter  connectivity, 
and  thus  topology,  and  since  topology  cannot  be  evaluated  locally,  there 
Is  a  strong  possibility  that  there  Is  no  low  level  test  which  can  be 
used  to  select  a  best  smoothing  mask.  In  general,  an  Ideal  smoothing 
mask  should  have  a  diameter  equal  to  the  width  of  the  line  whose 
connectivity  It  Is  enhancing.  If  the  mask  diameter  Is  smaller,  It  Is 
less  effective;  If  the  mask  diameter  Is  larger,  It  will  eventually 
eliminate  or  shift  the  line's  location  (see  Figure  12).  Ideally,  one 
should  separately  process  the  Image  with  each  of  the  available  smoothing 
masks  and  piece  together  the  best  delineations  after  a  final  evaluation. 
The  precise  control  structure  for  such  an  approach  Is  still  under 
Investigation.  The  alternative  we  are  using  at  present  Is  to  set  the 
mask  size  to  the  diameter  of  the  the  smallest  width  lines  we  wish  to 
ensure  are  Included  In  the  delineation.  In  our  experiments  we  used  the 
following  values: 
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D  ”  0  for  low  resolution  aerial  Images  where  line 
widths  of  1-2  pixels  are  common.  We  depend 
on  the  smoothing  already  Introduced  by  the 
limited  bandwidth  of  the  Image  acquisition  system. 

D  “  3  for  all  the  remaining  "normal"  Imagery. 

D  “  20  for  the  very  noisy  angiogram. 


e)  Comments:  For  a  given  line  of  nominally  constant  width,  the 
quality  of  Its  delineation  plotted  as  a  function  of  smoothing  mask 
diameter,  Is  a  parabolic  function  peaking  where  mask  size  equals  line 
diameter.  This  Is  the  only  method  we  have  found  so  far  to  get  a  good 
estimate  of  line  width. 


RIDGEPOINT  DETECTION 

a)  Purpose:  to  locate  a  set  of  points  which  contain  the  points 
comprising  the  linear  structures  In  an  Image.  These  points  should  be 
dense  along  the  rldge-llnes  (skeletons)  and  relatively  sparse  elsewhere. 
In  addition  to  location,  for  each  point  we  also  wish  to  obtain  a  measure 
of  "contrast"  between  the  point  and  the  background. 

b)  Algorithm:  rldgepolnts  (maxpolnts)  are  detected  by  searching 
for  local  (Intensity  or  other  attribute)  maxima  along  horizontal  and 
vertical  scan  lines.  Each  rldgepoint  Is  assigned  two  numerical  values 
which  are  measures  of  its  "local"  and  "global"  contrast.  These  values 
are  determined  as  shown  In  Figure  2,  and  the  larger  of  each  of  the  two 
values  determined  for  the  horizontal  and  vertical  directions,  are 
retained  for  each  of  the  two  attributes. 

c)  Parameter  Space:  with  the  possible  exception  of  scan-line 
positioning,  the  algorithm  for  rldgepoint  detection  has  no  Internal 
parameters. 

d)  Parameter  Selection:  A  dense  set  of  horizontal  and  vertical 
scan  lines  was  used  In  all  of  our  experiments.  Since  this  selection 
produced  excellent  results,  nothing  else  was  considered. 

e)  Comments:  The  competence  of  the  rldgepoint  algorithm  to 
abstract  the  linear  structure  of  an  Image  is  apparent  by  inspection  of 
Figures  4-10.  A  significant  feature  of  "ridepolnts”  is  that  their 
locations  are  Independent  of  any  monotonlc  Intensity  transformation  of 
the  Image,  and  their  geometric  configuration  Is  Independent  of  a  change 
of  scale. 
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ASSOCIATION  AND  THRESHOLDING 


a)  Purpose:  To  partition  the  ridgepoints,  detected  in  the  image, 
into  independent  coherent  subsets  at  some  maximum  level  of  complexity. 

b)  Algorithm:  Points  are  initially  grouped  into  mutually  exclusive 
clusters  such  that  every  point  in  a  cluster  is  within  some  maximum  2-D 
Eucledlan  distance  from  at  least  one  other  point  in  the  cluster.  The 
technique  used  is  a  fast  "two  pass"  algorithm  described  in  (Flrscheln 
and  Flschler  [5]).  If  the  number  of  points  in  a  cluster  is  greater  than 
some  preset  complexity  threshold,  then  the  cluster  is  thinned  by 
eliminating  those  points  with  the  lowest  contrast  (i.e.,  raising  the 
"contrast  threshold”);  a  separate  contrast  threshold  is  thus  determined 
Independently  for  each  cluster  in  a  image.  If  the  elimination  of  low 
contrast  points  causes  the  distance  condition  to  be  violated,  the 
cluster  is  decomposed  into  subclusters  which  now  assume  independent 
status,  and  the  above  procedure  is  iteratively  Invoked. 

c)  Parameter  Space:  The  256X256  image  contains  65,000  pixel 
locations,  and  the  total  number  of  ridgepoints  in  a  typical  scene  can 
vary  from  10,000  to  25,000  points.  An  initial  contrast  thresholding  is 
done  to  reduce  the  number  of  ridgepoints  to  about  10  percent  of  the 
image.  At  present,  the  global  contrast  threshold  is  always  set  to  be 
twice  the  value  selected  for  the  local  contrast  threshold.  Typical 
values  obtained  for  local  contrast  range  from  2  to  20  intensity  units 
for  a  picture  with  an  intensity  range  of  0  to  255  intensity  units. 

The  complexity  threshold  (the  maximum  number  of  points  in  a 
cluster)  is  based  on  both  the  competence  of  subsequent  techniques,  and 
on  computer  system  constraints.  It  is  currently  set  at  1500  points. 
The  contrast  thresholds  for  a  cluster  are  adjusted  by  the  algorithm  to 
satisfy  the  complexity  threshold  constraint. 

The  distance  threshold  employed  by  the  association  algorithm  is 
based  on  the  assumption  that  the  clearly  visible  linear  structures  in 
the  image  will  be  continuous,  and  only  a  slight  allowance  need  be  made 
for  acquisition  and  quantization  noise.  Values  of  2  and  3  are  the  only 
distances  considered;  this  allows  gaps  of  1  and  2  pixels,  respectively. 

d)  Parameter  Selection:  The  complexity  threshold,  with  a  setting 
of  1500,  allows  an3where  from  two  to  ten  times  more  ridgepoints  to 
appear  in  a  cluster  than  will  be  retained  in  the  final  delineation 
extracted  from  that  cluster;  it  is  extremely  conservative,  and  even 
doubling  or  halving  its  value  generally  has  almost  no  effect  on  the 
final  delineation.  The  distance  threshold  is  set  to  "2”  for  normal 
Images,  and  to  "3"  for  exceptionally  noisy  Images  which  have  been 
smoothed  with  a  very  large  diameter  filter  mask.  In  our  experiments,  a 
distance  threshold  of  2  was  used  for  all  the  images  except  the  angiogram 
where  we  used  a  distance  threshold  of  3. 
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LINKING 


a)  Purpose:  To  connect  the  elements  of  each  cluster  Into  a  graph 
(or  tree)  In  which  branches  delineate  potential  linear  structures  in  the 
image . 

b)  Algorithm:  Rldgepolnts  are  linked  (connected  to  their  nearest 
Eucledlan  distance  neighbors)  using  a  Minimum  Spanning  Tree  (MST) 
algorithm  (Relngold  [6]). 

c)  Parameter  Space:  There  are  no  Internal  parameters  In  the  MST 
algorithm. 


PRUNING 

a)  Purpose:  To  eliminate  rldgepolnts  which  do  not  fall  along  the 
major  linear  structures  In  the  Image. 

b)  Algorithm;  Based  on  the  assumption  that  Insignificant  linear 
segments  (and  accidental  alignments  of  noise  that  might  produce 
rldgepolnts)  will  result  In  short  terminal  branches,  or  produce  gaps 
near  the  terminal  ends  of  branches,  we  Iteratively  delete  from  the  tree 
terminal  branches  with  fewer  than  (1/3)X( ’’SHORT")  points,  and  those 
segments  at  the  ends  of  terminal  branches  which  have  a  gap  prior  to 
encountering  at  least  "SHORT”  contiguous  points;  when  no  additional 
points  can  be  deleted  by  further  Iterations,  those  terminal  branches 
with  fewer  than  "SHORT"  points  are  also  deleted. 

c)  Parameter  Space;  The  single  Internal  parameter  "SHORT"  Is  set 
to  the  length  of  the  shortest  linear  segment  we  would  be  willing  to 
consider  as  significant  in  the  final  output.  Without  knowing  anything 
about  the  Image  domain,  this  parameter  Is  based  on  purely  psychological 
considerations . 

d)  Parameter  Selection:  "SHORT"  has  been  hardwired  to  a  value  of 
15  for  256X256  Images;  this  value  gives  generally  good  results  both  In 
eliminating  noise  and  In  retaining  significant  structure. 

e)  Comments:  Eliminating  noise  segments  at  the  spanning  tree  level 
of  organization,  as  opposed  to  simply  eliminating  Isolated  or  low 
contrast  rldgepolnts,  Is  a  major  contributor  to  the  excellent 
performance  of  the  complete  delineation  algorithm.  In  an  experiment 
with  20  randomly  generated  binary  overlays,  representing  detected 
rldgepolnts  (rldgepolnt  density  lOZ),  no  randomly  formed  linear  segments 
survived  this  filtering  step. 
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LINEAR  SEGMENT  EXTRACTION 


a)  Purpose;  To  extract  Individual  coherent  linear  segments  from 
the  MST. 

b)  Algorithm:  It  Is  assumed  that  the  MST  for  a  given  cluster, 
after  the  pruning  step,  may  contain  a  number  of  Inter-twlned 
perceptually  significant  linear  segments.  The  problem  Is  to  determine 
at  each  node  of  the  MST  which  branches  should  be  linked  to  each  other, 
and  which  branches  should  be  split  off  as  end-point  (terminal)  portions 
of  separate  linear  segments.  The  algorithm  considers  each  node  of  the 
tree  In  some  arbitrary  order;  at  each  node  It  measures  four  attributes 
of  each  Incoming  pair  of  branches:  the  "angle"  the  branches  form  with 
their  vertex  at  the  given  node,  the  difference  between  the  average 
intensity  grayscale  values  of  the  branches,  whether  both  branches  are 
terminal  branches,  and  whether  either  of  the  branches  has  a  gap 
Immediately  adjacent  to  the  given  node.  A  "mismatch"  score  is  computed 
for  each  combinatorially  selected  pair  of  branches,  and  that  pair  with 
the  lowest  score  is  designated  the  "control  pair." 

If  the  score  computed  for  the  control  pair  is  greater  than  a 
locally  computed  "compatibility  threshold,"  the  node  will  be  deleted 
from  the  MST,  and  the  MST  thus  decomposed  Into  a  number  of  smaller 
disjoint  trees;  none  of  the  branches  entering  the  given  node  will  be 
connected  to  each  other  in  the  final  delineation. 

If  the  mismatch  score  computed  for  the  control  pair  is  less  than 
the  compatibility  threshold,  then  the  control  pair,  and  any  other 
segment  which  when  paired  with  one  of  the  segments  in  the  control  pair 
has  a  score  almost  as  good  as  that  of  the  control  pair  (within  5%)  will 
remain  connected  to  the  given  node.  All  other  branches  will  be  split 
off  from  the  given  node  but  will  be  joined  at  a  corresponding  node  of  a 
new  tree  and  then  subjected  to  the  procedure  just  described  (this  allows 
more  than  one  pair  of  branches,  entering  a  given  node,  to  be  linked  into 
a  continuous  path). 

After  all  of  the  nodes  have  been  subjected  to  the  above  procedure, 
the  maximum  length  paths  in  all  of  the  resulting  trees  will  be  extracted 
as  independent  linear  segments.  Branches  split  off  from  the  maximum 
length  paths  form  new  trees  which  are  recursively  processed  in  the  same 
way. 


Parameter  Space:  This  procedure  has  five  parameters;  these  are  the 
compatibility  threshold  and  weighting  factors  assigned  to  each  of  the 
attributes  in  computing  the  incompatibility  score  for  a  pair  of  branches 
entering  a  given  node.  All  attributes  are  computed  using,  at  most,  the 
eight  points  closest  to  the  node  along  each  of  the  entering  branches. 

a)  Compatibility  Threshold:  computed  for  each  node  as  20%  of  the 
average  grayscale  value  (AGSV)  of  all  the  branches  entering  the  node. 
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b)  Penalty  for  two  branches  meeting  at  an  acute 
angle  =  (6%)X(AGSV). 

c)  Reward  for  two  branches  meeting  at  an  angle  of  at  least  145 
degrees  =  (8%)X(AGSV). 

d)  Penalty  for  both  branches  being  terminal 

branches  =  (5%)X(AGSV). 

e)  Penalty  for  either  (or  both)  branches  having  a  gap  immediately 
adjacent  to  the  given  node  “  (4%)X(AGSV). 

f)  Penalty  for  difference  in  average  gray  scale  value  of  the  two 
branches  =  (actual  gray  scale  difference). 

Parameter  Selection:  The  problem  of  measuring  local  similarity,  and 
best  geometric  continuation,  of  a  set  of  line  segments  entering  a  common 
node,  is  at  a  level  of  complexity  requiring  a  rule-based  system  to  make 
a  determination  in  reasonable  agreement  with  human  perceptual 
judgements.  The  particular  attributes  chosen  were  ones  which  seemed 
relevant  and  could  be  easily  measured  (line-width  would  have  been 
Included  as  an  attribute  if  we  had  a  good  way  of  measuring  it  locally; 
our  method  of  measuring  the  "angle”  between  two  segments  will  be 
improved  in  the  near  future).  The  selected  values  of  the  weighting 
factors  were  based  on  informal  experiments,  and  then  hard-wired  into  the 
system. 

Comments:  It  is  obvious  that  the  selected,  attributes,  our  method  of 
measuring  them,  and  the  assigned  weighting  factors,  can  all  be  Improved 
by  a  more  systematic  study  of  their  psychological  implications; 
nevertheless,  the  currently  Implemented  procedure  seems  to  make  very  few 
poor  decisions  and  does  not  appear  subject  to  major  improvement. 


LINEAR  SEGMENT  DECOMPOSITION 

Purpose;  To  split  individual  linear  segments  that  are  not 
perceptually  "homogenous." 

Algorithm;  The  preceding  procedure  for  linear  segment  extraction 
accomplishes  this  operation  when  two  or  more  segments  come  together  at  a 
common  point  (node).  However,  it  is  possible  that  exactly  two  distinct 
segments  can  accidentally  meet  at  their  nominal  end-points,  and  a 
special  procedure  may  be  required  to  check  for  this  situation.  No 
algorithm  is  currently  included  for  this  purpose.  However,  the  linear 
segment  extraction  algorithm  could  be  applied  at  every  point  along  a 
given  line  segment  to  achieve  this  purpose  if  desired. 
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LINEAR  SEGMENT  EVALUATION 


Purpose;  To  rank  order  a  given  set  of  linear  segments  with  respect 
to  perceptual  quality  (e.g.,  given  two  alternative  potential 
delineations,  which  one  Is  more  acceptable  from  a  perceptual  standpoint; 
or  assuming  we  wish  to  restrict  our  "sketch"  of  an  Image  to  the  "best"  N 
linear  segments,  how  the  selection  Is  to  be  made). 

Algorithm:  This  procedure  separately  measures  two  attributes  of 
each  line  segment  —  Its  photometric  quality  (contrast  with  the 
background)  and  Its  geometric  quality  (coherent  length  and  smoothness). 
Photometric  quality  Is  measured  by  simply  summing  the  Local  Maxima 
values  assigned  to  each  point  of  the  segment.  Geometric  quality  Is 
computed  by  adding  the  number  of  points  In  the  segment  to  the  distance 
between  end  points  of  the  segment  and  then  subtracting  three  times  the 
number  of  gaps,  and  also  the  number  of  locations  at  which  the  segment 
has  discontinuity  In  local  smoothness.  Numerical  values  are  normalized 
by  converting  them  Into  rank  orderings  for  all  of  the  segments  being 
evaluated.  The  final  score  assigned  to  each  segment  Is  the  sum  of  Its 
photometric  and  geometric  orderings,  with  the  geometric  score  weighted 
twice  that  of  the  photometric  score. 

Parameter  Space:  Four  quantities  are  measured  for  each  line 
segment;  contrast,  number  of  points,  number  of  gaps,  and  number  of 
discontinuities  In  smoothness.  The  final  ranking  Is  a  weighted 
combination  of  the  values  of  these  attributes. 

Parameter  Selection:  The  relative  weightings  assigned  to  the 
attributes,  as  described  above,  were  determined  by  Informal 
experimentation  and  the  hard-wired  Into  the  system. 

Comments:  More  than  any  other  of  the  preceding  steps  In  the 
delineation  process,  evaluation  must  be  based  on  psychological  factors, 
and  thus  effective  performance  depends  on  observing  and  modeling  human 
behavior  In  this  task.  Our  presently  Implemented  algorithm,  described 
above.  Is  a  first  unsophisticated  attempt  to  achieve  this  goal  --  Its 
main  defects  are  lack  of  an  effective  way  to  measure  local  width  and 
local  smoothness.  This  particular  procedure,  while  not  currently 
essential  to  the  overall  performance  of  the  delineation  algorithm,  can 
play  a  vital  role  In  automatically  adjusting  the  heurlstlcally  set 
parameters  of  the  other  procedures  through  performance  evaluation  feed¬ 
back.  In  order  to  take  advantage  of  such  an  approach,  the  psychological 
model  underlying  this  procedure,  and  the  methods  for  measuring  the 
relevant  attributes,  still  need  significant  Improvement.  Other  work  now 
underway  Is  addressing  some  of  these  Issues  (reference  Flschler  and 
Bolles  I?]). 
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{c!  SMOOTHED  GRAY  SCALE  IMAGE  (d)  INTENSITY  PROFILES  OF  SMOOTHED  IMAGE 


(e)  terrain  map  OF  INTENSITY  SURFACE  OF  (f)  LINEAR  FEATURE  POINTS  FOUND  IN  VALLEYS 

SMOOTHED  IMAGE  OF  INTENSITY  SURFACE 

FIGURE  1  IMAGE  SMOOTHING,  AND  FINDING  LINEAR  FEATURE  POINTS  IN  THE  INTENSITY 
VALLEYS  OF  A  RADIOGRAPHIC  IMAGE  (ANGIOGRAM) 


(b)  ORIGINAL  GRAY  SCALE  IMAGE 


(b)  INTENSITY  PROFILES  OF  ORIGINAL  IMAGE 
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LLM  =  left  local  minimum  Intensity 

LGM  =  left  global  minimum  intensity 

RLM  =  right  local  minimum  intensity 

RGM  =  right  global  minimum  intensity 

LIM  =  "local  maxima"  of  intensity  difference 

GIM  =  "global  maxima"  of  intensity  difference 


Point 

Point  intensity 

LLM 

RLM 

LGM 

RGM 

LIM 

GIM 


FIGURE  2  EXAMPLES  OF  LOCAL  AND  GLOBAL  MAXIMA  COMPUTATION 
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Q 


(a)  A  SINGLE  CLUSTER  OF  LINEAR  FEATURE  POINTS 


(b)  A  SINGLE  CLUSTER  OF  LINEAR  FEATURE  POINTS 


(e)  THE  TWO  LINEAR  SEGMENTS  EXTRACTED  (f)  THE  LINEAR  SEGMENT  EXTRACTED 

{one  horizontsl  line  and  one  vertical  line) 

FIGURE  3  PRUNING  AND  LINEAR  SEGMENT  EXTRACTION  FROM  A  SINGLE  CLUSTER 
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(a)  ORIGINAL  GRAY  SCALE  IMAGE:  DROAD2  DETECTED  LINEAR  FEATURE  {INTENSITY 

MAXIMA)  POINTS 


{c)  THRESHOLDED  LINEAR  FEATURE  POINTS  W)  SEGMENTS  FOUND  IN  THRESHOLDEO  LINEAR 

FEATURE  POINTS 


(e)  SEGMENTS  WITH  30  OR  MORE  POINTS 


(f)  10  HIGHEST  RANKED  SEGMENTS 


FIGURE  4  EXAMPLE  SHOWING  THE  SEQUENCE  OF  STEPS  IN  THE  DELINEATION 
PROCESS  FOR  A  LOW  RESOLUTION  AERIAL  IMAGE  (DROAD2j 
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(a)  ORIGINAL  GRAY  SCALE  IMAGE;  AFROAD1 


<b)  DETECTED  LINEAR  FEATURE 
(INTENSITY  MAXIMA)  POINTS 


(e)  SEGMENTS  WITH  30  OR  MORE  POINTS  (t)  9  HIGHEST  RANKED  SEGMENTS 


FIGURE  5  EXAMPLE  SHOWING  THE  SEQUENCE  OF  STEPS  IN  THE  DELINEATION 
PROCESS  FOR  A  LOW  RESOLUTION  AERIAL  IMAGE  (AFROAD1) 
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(c)  THRESHOLDEO  LINEAR  FEATURE  POINTS  SEGMENTS  FOUND  IN  THRESHOLDED 

LINEAR  FEATURE  POINTS 


{c)  SEGMENTS  WITH  30  OR  MORE  POINTS  (f)  27  HIGHEST  RANKED  SEGMENTS 


FIGURE  6  EXAMPLE  SHOWING  THE  SEQUENCE  OF  STEPS  IN  THE  DELINEATION 
PROCESS  FOR  A  HIGH  RESOLUTION  AERIAL  IMAGE  (FARM-ULl 
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(c)  THRESHOLDEO  LINEAR  FEATURE  POINTS  W)  SEGMENTS  FOUND  IN  THRESHOLOED 

linear  feature  points 


{a!  SEGMENTS  WITH  30  OR  MORE  POINTS  (f)  10  HIGHEST  RANKED  SEGMENTS 


FIGURE  7  EXAMPLE  SHOWING  THE  SEQUENCE  OF  STEPS  IN  THE  DELINEATION 
PROCESS  FOR  A  HIGH  RESOLUTION  AERIAL  IMAGE  {MARtN2) 
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(a)  ORIGINAL  GRAY  SCALE  IMAGE;  PCB1 


(b)  DETECTED  LINEAR  FEATURE 
(INTENSITY  MAXIMA)  POINTS 


(cl  THRESHOLDED  LINEAR  FEATURE  POINTS 


(e)  SEGMENTS  WITH  30  OR  MORE  POINTS 


(d)  SEGMENTS  FOUND  IN  THRESHOLDED 
LINEAR  FEATURE  POINTS 


(f)  24  HIGHEST  RANKED  SEGMENTS 


FIGURE  8  EXAMPLE  SHOWING  THE  SEQUENCE  OF  STEPS  IN  THE  DELINEATION 
PROCESS  FOR  A  PRINTED  CIRCUIT  BOARD  IMAGE  (PCBI) 
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(a)  ORfGINAL  GRAY  SCALE  IMAGE:  ANG105A 


(b)  DETECTED  LINEAR  FEATURE 
(INTENSITY  MAXIMA)  POINTS 


(e)  SEGMENTS  WITH  30  OR  MORE  POINTS 


(f)  6  HIGHEST  RANKED  SEGMENTS 


FIGURE  9  EXAMPLE  SHOWING  THE  SEQUENCE  OF  STEPS  IN  THE  DELINEATION 
PROCESS  FOR  A  RADIOGRAPHIC  ANGIOGRAM  (ANGI05A) 


(al  ORIGINAL  GRAY  SCALE  IMAGE:  TURTLE  (b)  EXTRACTED  LINE  SEGMENT  FOR  TURTLE 


{a}  ORIGINAL  GRAY  SCALE  IMAGE:  TREE  (d)  EXTRACTED  LINE  SEGMENT  FOR  TREE 


(e)  ORIGINAL  GRAY  SCALE.  !R?AGE:  MOUSE  (f)  EXTRACTED  LINE  SEGMENT  FOR  HOUSE 


FIGURE  TO  EXTRACTION  OF  LINEAR  SEGMENTS  FROM  PENCIL  DRAWINGS 
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a.  Gray  scale  Intensities  of  a  thick,  noisy,  horizontal  line  In  window 
237,169  253,183  of  unsmoothed  Image,  PCBl. 

Line  has  Intensity  values  greater  than  200. 

*-ed  pixels  Indicate  points  at  which  global  maxima  Is  greater  than  15 


237  238  239  240  241  242  243  244  245  246  247  248  249  250  251  252  253 

169  164  168  172  175  177  178  179  179  179  180  180  179  177  174  170  167  166 

170  176  177  179  180  181  181  181  181  181  181  181  181  180  178  176  175  175 

171  195  195  195  195  195  195  195  195  194  194  194  193  193  191  190  190  190 

172  212  212  212  211  211  211  211  211  211  210  210  209  208  207  206  206  206 

173  223  223  223  223  223  223  223  223  223  222  221  221  220  219  218  218  218 
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b.  Gray  scale  Intensities  after  smoothing  with  mask  of  diameter  11. 

*-ed  pixels  Indicate  points  at  which  global  maxima  Is  greater  than  15 


FIGURE  11  SMOOTHING  WITH  AN  APPROPRIATE  MASK  WILL  ENHANCE 
LINE  STRUCTURE 
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a.  Gray  scale  Intensities  of  parts  of  two  segments  in  window  167,43 
183,57  of  unsmoothed  image,  DR0AD2. 

*-ed  pixels  Indicate  points  on  two  distinct  linear  segments. 


167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 

180 

181 

182 

183 

43 

110 

112 

113*114 

113 

111 

107 

103 

100 

97 

94 

91 

87 

83 

79 

76 

74 

44 

115 

118 

120*122 

121 

120 

117 

113 

110 

107 

103 

99 

94 

89 

84 

80 

76 

45 

118 

122 

126 

128*129 

129 

127 

124 

121 

118 

114 

109 

103 

97 

90 

85 

80 

46 

118 

124 

129 

133 

135*136 

135 

133 

131 

128 

125 

119 

113 

106 

98 

92 

86 

47 

117 

123 

129 

134 

138 

140*140 

140 

139 

136 

133 

128 

122 

115 

107 

99 

92 

48 

114 

120 

126 

133 

138 

141*143*143 

143 

141 

138 

134 

128 

122 

115 

107 

99 

49 

113 

118 

124 

130 

136 

140 

143*144*145 

143 

141 

138 

133 

128 

121 

114 

106 

50 

112 

116 

122 

128 

134 

139 

142 

144*145 

144 

143 

140 

137 

132 

126 

120 

111 

51 

113 

116 

121 

126 

132 

137 

140 

143*145*145 

144 

142 

139 

135 

130 

124 

116 

52 

114 

116 

120 

125 

130 

135 

139 

142 

143*144 

144 

142 

140 

137 

132 

127 

119 

53 

115 

117 

120 

124 

129 

133 

137 

140 

142*143 

143 

142 

141 

138 

134 

129 

121 

54 

118 

118 

120 

123 

128 

131 

135 

137 

140 

141*142 

142 

141 

139 

135 

131 

124 

55 

121 

120 

121 

123 

126 

130 

132 

135 

137 

139 

140*140 

140 

139 

136 

132 

125 

56 

124 

122 

122 

123 

125 

127 

130 

132 

134 

136 

138 

139*139 

139 

137 

133 

127 

57 

128 

124 

122 

123 

124 

125 

127 

128 

131 

133 

136 

137*138 

138 

137 

134 

129 

b.  Gray  scale  intensities  after  smoothing  with  mask  of  diameter  21. 
*-ed  pixels  indicate  points  on  the  linear  segment  found. 

Upper  right  segment  has  been  eliminated. 


FIGURE  12  SMOOTHING  WITH  A  LARGE  DIAMETER  MASK  MAY  ELIMINATE  OR 
SHIFT  SEGMENTS 
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